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Abstract. We calculate the squeezed limit of the bispectrum produced by inflation with 
multiple light fields. To achieve this we allow for different horizon exit times for each mode 
and calculate the intrinsic field-space three-point function in the squeezed limit using soft- 
limit techniques. We then use the 6N formalism from the time the last mode exits the 
horizon to calculate the bispectrum of the primordial curvature perturbation. We apply our 
results to calculate the spectral index of the halo bias, nsb, an important observational probe 
of the squeezed limit of the primordial bispectrum and compare our results with previous 
formulae. We give an example of a curvaton model with ngb ~ 0{ns — 1) for which we 
find a 20% correction to observable parameters for squeezings relevant to future experiments. 
For completeness, we also calculate the squeezed limit of three-point correlation functions 
involving gravitons for multiple field models. 
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1 Introduction 

Reliable calculations of the N-point correlation functions of the primordial curvature per¬ 
turbation, are essential to confront models of inflation with present and future observa¬ 
tional constraints. Since inflation can occur at energies as high as lO^^GeV, these correlation 
functions provide an unparalleled observational window into high energy physics, providing 
information about the helds and their interactions active in the early universe. In this paper 
we consider the three-point correlation function, and discuss how it can be calculated in a 
particular limit, known as the ‘squeezed limit’. Our primary aim is to provide clarity on how 
to accurately confront models of inflation with more than one light held against observations 
sensitive to this limit. This includes models in which more than one held supports inhation. 
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as well as spectator models such as the curvaton scenario and modulated reheating models 
- we refer to such models as multi-field models, and exclude cases in which additional heavy 
fields play a role. 

The three-point function, parametrised by the bispectrum, is a function of three wave 
vectors which sum to zero as a result of momentum conservation, forming a triangle in momen¬ 
tum space. The squeezed limit refers to the case where one of the associated wave numbers is 
much smaller than the other two, such that the triangle looks ‘squeezed’. Calculations of the 
three-point function are now extremely mature, yet for technical reasons previous multiple 
field calculations have only been been performed explicitly for the case of a mild hierarchy 
between wave numbers, as we will see. The highly squeezed limit is, however, very important 
both from an observational and a theoretical point of view. 

The squeezed limit of the the bispectrum is the simplest possible example of a more 
general class of limits of correlation functions, referred to as soft limits. Soft limits occur 
when there exists a separation of scales in a physical problem. In the inflationary context, 
soft limits of correlation functions of the primordial curvature perturbation offer an exciting 
opportunity to confront theory with observations. For example, in the case of a single slow-roll 
field with canonical kinetic terms, Maldacena [1] found that the bispectrum of the curvature 
perturbation in the squeezed limit is purely determined by the tilt of the power spectrum, 
with the assumption of a Bunch-Davies initial state. The relation is: 12/^^ = —^i^s — 1), 
where is the squeezed limit of the reduced bispectrum, and — 1 = —0.032 ± 0.006 [2] 
is the spectral tilt. Creminelli & Zaldarriaga [3] (see also Ref. [4]), showed that Maldacena’s 
result holds even without the assumption of slow-roll^ in all models with a Bunch-Davies 
initial state and where the classical solution is a dynamical attractor - the proof of which 
was later formalized by Cheung et al. [7]. Thus a detection of /nl~ 0(0.01) would rule out 
all single field models with a Bunch-Davies initial state and where the classical solution is a 
dynamical attractor^. For single field inflation, considerable work has also gone into studying 
more general soft limit results and providing further consistency relations amongst correlation 
functions [8, 10-31]. Moreover, soft limits can be used to provide information about other 
fields present during inflation [32-39]. 

In the case of inflation with multiple light fields, in contrast to the single field case, model 
independent results such as the Maldacena consistency relation are not possible. However, 
many observables which constrain non-Gaussianity produced by multiple field inflation are 
particularly sensitive to soft limits. Examples include^ the spectral index of the halo bias [45] 
and CMB //-distortions [46]. This is why explicit calculations of soft limits in the multiple 
field context are important when comparing model dependent predictions against observation. 

In this paper, as a starting point to more general studies of soft limits in multiple field 
inflation, we explore the simplest case of the squeezed limit of the bispectrum. In most 
previous studies, the path to the bispectrum for multiple field models has been to first use 
the in-in formalism to calculate the three-point function of scalar field perturbations at a time 
soon after all modes have left the horizon, as was first done by Seery &: Lidsey [47]. Next 
the 5N formalism [48, 49] is applied to convert the field-space correlations to correlations 

^See, for example, [5, 6] for single field, but non-slow-roll models which obey Maldacena’s relation. 

■^See [8] (and e.g. [9]) where more general initial states are considered, and see [10] (and references therein) 
where non-attractor models are considered. 

^One might naively think that the long-wavelength fluctuations used to model the observed CMB power 
asymmetry (see e.g. [40-44]) may also be an observational probe of soft limits - however, to describe the 
asymmetry, the long wavelength mode is required to be superhorizon, and so won’t be suitable as a soft 
momentum in the correlation functions considered in this paper. 
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of ( [50]. The result of Seery &: Lidsey for the three-point function, however, requires that 
there is not a large hierarchy between the three wavenumbers involved in the bispectrum, 
and thus that the modes of the bispectrum must cross the horizon during inflation at roughly 
the same time. As is clearly stated in their paper, therefore, their result is not valid in the 
highly squeezed limit where there is an appreciable difference in the exit times of different 
modes. Moreover, on using SN to convert from held-space fluctuations to C, one hnds that the 
three-point function of the curvature perturbation involves copies of the two-point function 
of held huctuations evaluated after all modes have exited the horizon. At this point in the 
procedure, previous explicit calculations have considered at most a mild hierarchy between 
the scales at which these two-point functions are evaluated [51, 52], with | log(A;i/A: 3 )| 0(1), 

where ki is the long-wavelength mode, and is the short-wavelength mode. 

Future experiments, together with the expected amount of squeezing they will be sen¬ 
sitive to, are shown in Table 1 [52]. These experiments will probe a hierarchy much larger 
than that allowed by previous theoretical calculations. Therefore, it is very important to have 
theoretical predictions for multiple held models valid in the highly squeezed limit, to be able 
to compare with observations. 


Experiment 

Dark Energy Survey 

Euclid 

/i-distortions 

Squeezing 

log(/ci//c3) ~ -2 

log(fci/fc3) ~ -8 

log(fci/fc3) ~ -19 


Table 1: Future experiments and their observable range of scales. 


In our study, therefore, we wish to relax the requirement of a mild hierarchy, and study 
models explicitly in the highly squeezed limit. To do so we will use a similar soft limit 
argument to Cheung et al. [7], but applied to calculate the three-point function of the scalar 
held perturbations. This result reduces to that of Seery & Lidsey for mild squeezing, but is 
valid in the highly squeezed limit and does not rely on slow-roll. In analogy with the single 
held case, where the three-point function depends on the tilt of the power spectrum, we hnd 
that the three-point function of the scalar held perturbations depends on derivatives of the 
held space two-point function with respect to the background value of the helds. Armed with 
this result, we further relax the usual assumption that the copies of the two-point function 
of the held perturbations, which appear in the three-point function for ([, involve only a mild 
hierarchy of scales. This can easily be achieved by accounting for the evolution between 
horizon crossing times. We hnd that for any model in which previous applications of the 
5N formalism give rise to analytic results, analytic expressions for the highly squeezed limit 
are also possible. We hnish by considering a specihc curvaton model and compare our new 
squeezed limit formulae with previous expressions. We hnd signihcant differences for cases in 
which non-Gaussianity depends on scale. 

The outline of this paper is as follows: in §2 we review the previous applications of the 5N 
formalism, taking a pedagogical approach in this section which will be helpful when we come 
to extend previous work. In §3 we conveniently parametrise the evolution of the superhorizon 
held perturbations between exit times in terms of a ‘F-matrix’, and calculate the three-point 
function of the scalar held perturbations at the time the last mode exits the horizon using 
soft-limit arguments. Putting these elements into the 5N formalism we calculate the highly 
squeezed limit of the bispectrum of the curvature perturbation. In §4 we calculate the scale 
dependence of the squeezed limit of the reduced bispectrum, focussing on the spectral index 
of the halo bias. In §5 we provide explicit formulae for the F matrix, and investigate the 
concrete example of the mixed inhaton-curvaton scenario [53] with self-interactions [43]. 
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Throughout this paper we work in units where h = c = 1 and we set the reduced Planck 
mass Mp = 1. 

2 Review oi 5N 

In this section we review previous applications of the 5N formalism to calculate the bispectrum 
of Q. In doing so, we present key dehnitions and formulae essential for the subsequent sections. 

2.1 Defining C, 

We will be considering Fourier modes which make up a given triangle of the bispectrum, 
labelled by three wave vectors ki, \^2 and ks. These scales satisfy the condition ki+k 2 +k 3 = 
0 because of momentum conservation, and without loss of generality we order them such that 
ki < k 2 < ks, where k = |k|. We denote the times these modes cross the horizon by {ti,t 2 ,t 3 }, 
dehned by the condition ki = a{ti)H{ti) etc. and so we have ti <t 2 < For the bispectrum 
there are only two limiting possibilies for the relative magnitudes of wave numbers. The close 
to the equilateral case, with fei ~ A :2 ~ fea, and the squeezed shape, with ki <C A :2 ~ /cs. The 
amount of squeezing can be measured by the value of log(/ci//c 3 ), a negative number with 
greater absolute value indicating greater squeezing. 

In the context of multiple held inhation, the cosmological perturbation of most interest 
is the uniform-density curvature perturbation, which coincides with the comoving curvature 
perturbation on super horizon scales, where it is conserved if the system reaches an adiabatic 
limit [54, 55]. In general the curvature perturbation can be dehned as a scalar perturbation 
to the spatial metric for a given foliation of spacetime, and moreover can be written as a local 
perturbation to the scale factor as [54, 56, 57] 

a(t,x) = (2.1.1) 

Choosing the t-slicing to be such that the spatial hyper surfaces have uniform-density (UD) 
leads to V'ud(Tx) = while a hat t-slice is dehned by V’flat(Cx) = 0- The power 

spectrum, P^, and bispectrum, of the curvature perturbation are then dehned by the 
two and three-point correlation functions 

(CkiCk.) = Pdki){27rf6{ki + ka) (2.1.2) 

(CkiCk2Ck3) = Bc;{ki, k2, /C3)(27r)^h(ki +k2 + ks). (2.1.3) 

2.2 The 6N formalism 

Using the spatially dependent dehnition of the scale factor, the number of e-folds which occurs 
between two time slices of the perturbed spacetime, labeled by T and tu respectively, is a 
function of position and is given by 

Pu nf-f Pu 

N{T,tu,^)^ / .,’ ^ dt = / R(t)dt-h V'(i«,x) - V'(T,x) (2.2.1) 

Jt o(t,x) Jt 

while the unperturbed number of e-folds is given by iVo(T, t„) = H{t)dt. Taking the 
T-slices to be hat, and the t^-slices to be uniform density gives 

C(4,x) = N{T,tu,^) - No{T,tu) = (5iV(t«,x) 
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( 2 . 2 . 2 ) 



which is the ceiebrated 6N formuia [48]. We note that 5N doesn’t depend on the initiai 
time T [54], Typicaiiy T is taken to be some time after aii the modes invoived in a given 
correiation function have exited the horizon (for the bispectrum this means T > t^). If the 
system becomes adiabatic C is also independent of tu and so the later time should be the 
time at which adiabaticity is reached. If it doesn’t become adiabatic, then tu can be taken to 
be the time at which the correlations are required (see for example the discussion in Ref. [58]). 

To use this formalism in practice, we must employ the separate universe approximation 
[49, 59] to cosmological perturbation theory. This states that on super-horizon scales positions 
in the perturbed universe evolve independently of one another, and do so according to the 
same equations as the unperturbed cosmology, so that every position can be treated as a 
‘separate universe’. The number of e-folds which occur at every position can therefore be 
calculated using the local conservation and Friedmann equations in that separate universe. 

For an inflationary model with n scalar fields, where i runs from one to re, we 
can split the held values on any flat slice into background and perturbed parts = 

^i{t) + 5(j)i{t,x.). Further demanding that the slow-roll equations of motion are satisfied at 
time T, such that 3H^ = V{(j)i) and SFfc/ij = —Vi , the initial conditions at time T for the 
perturbed cosmology become dependent only on (j)i (T, x). Consequently, the number of e- 
folds to any subsequent time slice becomes a function only of the initial held values, even if 
the cosmology evolves away from slow-roll. In particular, we can write 

N{T,tu,^) = N{{^,{T,^)),tu) = N{{MT) + 6MT,^)),tu) (2.2.3) 

which gives 

C{x) = N{{cPi{T) + 6(t>i{T,^)),tu) - No{UT),tu). (2.2.4) 


One can then Taylor expand in the initial flat slicing field perturbations, which in Fourier 
space leads to [48, 50] 


Cki 

where 



(2.2.5) 

( 2 . 2 . 6 ) 


and where * denotes convolution. The additional vector (boldface) subscript indicates wavevec- 
tor and the superscript in brackets is a shorthand indicating evaluation time 

= MT) (2.2.7) 


which we use from now on. 

The correlation functions of the held perturbations are defined as 

= sjf (A:i)(27r)3h(ki +k2) 

iki,k2, fc3)(27r)3<5(ki + k 2 + kg) 

such that the power spectrum [48] and bispectrum [50] of () are given by 

Pc(fca) =NPNp^^\ka) 

h) h) 

+ ivf ) [sjf (A:i)s(P(fc 2 ) + (fei ^k 2 ^ k,)] 


( 2 . 2 . 8 ) 

(2.2.9) 


( 2 . 2 . 10 ) 

( 2 . 2 . 11 ) 
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where the arrows indicate there are two additional terms formed by cyclic permutations. 

One can then dehne the reduced bispectrum, conventionally denoted by /NL(fci) ^ 2 ,^: 3 ), 
by comparing the amplitude of the bispectrum to the square of the power spectrum as follows 


fNL{h,k2,k3) 


5 Bc_{ki,k2,k3) __ 

6 [Pt;{ki)P(^{k2) + (fci ^2 ks)]' 


( 2 . 2 . 12 ) 


2.3 5N for /ci ~ /c 2 ~ k^ 

For light, canonically normalized helds, the held-space correlation function for a given wavenum¬ 
ber takes a very simple form at the time the wavenumber crosses the horizon^ (and from which 
the 5N formalism can be employed), one hnds [60, 61] 

( 1 ) 

4 (*■■>) = 

In the most common application of 6N, one assumes a near-equilateral momentum regime, 
where all the wavenumbers are approximately equal, fei ~ A ;2 ~ k^ and thus the horizon 
crossing times of the three wavenumbers involved in the bispectrum can be identihed with a 
single time, t*, such that ss t 2 ~ ^3 ~ In this regime, it is then common to make the 
simple choice T = t*. 

For canonical slow-roll inflation, one can take^ 

( 2 . 3 . 2 ) 


and similarly for the other wavenumbers, k 2 and k^. Moreover, Seery & Lidsey [47] used the 
in-in formalism to calculate Oi^ijl.{ki, A: 2 , ks) in the close to equilateral momentum conhguration 
at the time t*. They found 



{ki,k2,k3) 


/^O 1/^0 

^1^2^3 


^ Hjk 
27 r I ^ 4 i 7 (*) 

/ 6 perms 


(_ 3 MM 


2.2^2 

k? 


{ki + 2 /^ 3 ) -|- —kf — 



( 2 . 3 . 3 ) 


where the sum is over the six permutations of (ijk) while simultaneously rearranging the 
momenta ki,k 2 ,k 3 such that the relative positioning of the fc’s is respected. We emphasize 
that this result assumed that the three k modes crossed the horizon at roughly the same time 
and cannot be trusted when the crossing times are too different®, as in the case of the highly 
squeezed limit, which we will consider in this paper. 

Strictly speaking this result is the form the two-point function takes once the decaying mode present at 
horizon crossing has been lost, written in terms of horizon crossing parameters. 

®The leading order in slow-roll correction to this [60-62] is to replace 5ij 1 —>■ 5ij +2cuij , where c = 2—log 2 — 7 , 
with 7 the Euler-Masheroni constant, and Uij = —(logF),ij, but we will not use this correction in this work. 

®This is because in evaluating the time integrals of the in-in calculation, the time-dependent coefficients 
of the field perturbations were all evaluated at the common horizon crossing time, t*, of the near-equilateral 
modes. This can’t be done if the modes exit at largely different times, as is the case in the highly squeezed 
limit. 
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Sticking with the near-equilateral conhguration, with fci ~ A ;2 ~ and taking the 6N 
initial time T = the results (2.3.1) and (2.3.3) can then be used together with Eqs. (2.2.10)- 
(2.2.12) to give the well-known expression [50], valid when /NL(fc*) ^(f)) 


/nl 


^ J 

6 (lv/*)iv/*^)2 


(2.3.4) 


where (2.3.3) has been used to demonstrate that the second term in Eq. (2.2.11) must domi¬ 
nate over the hrst term if /nl is to be large - i.e. the contribution from a is neglected [63]. 
Eq. (2.3.4) retains a dependence on fc* through the horizon crossing time t*, and the relation 
that 

Given that it retains this dependence on the crossing scale, which is a measure of the 
overall scale of the bispectrum in the near-equilateral conhguration, the bispectrum calcu¬ 
lated above is sometimes referred to as quasi-local (to be contrasted with the local shape, in 
which /nl is independent of all three fc’s). Differentiating (2.3.4), Byrnes et al. [51] gave an 
expression for the tilt of /nL) denoted equilateral triangles 


n> = 

/NL 


where n) ' = 


(*) _ dlog|/NL| _ 


dlog A: 

1 dlogA^f^ 

H dt* 


(*) , (*) 
+ ny 


' 27 / ^ 27 7 




— 4h 


(*) _ 

= 


dlogAT^^^ 


H dt* 


and h = 


jyi*) y*) 
J J 


(2.3.5) 

(2.3.6) 


2.4 Beyond the near-equilateral configurations 


In the present paper we aim to be go beyond the case where A;i ~ A :2 ~ A: 3 , to study the 
highly squeezed limit. Before we do, we must stress that although most previous studies use 
the formulae presented in the previous sub-section, some attempts have already been made to 
relax this near-equilateral assumption, which we briehy review now. Full details will appear 
when we compare our expressions to earlier ones in later sections. 

In Ref. [51], in addition to the consideration of scale dependence discussed above, Byrnes 
et al. considered the more general dependence of the reduced bispectrum on ki, k 2 and k^, 
using the 5N framework. They did so by solving the slow-roll held equations for scalar 
huctuations to hrst order in the quantity | log(A:i/A: 3 )| about the time t*. This was used to 
allow a mild hierarchy of scales. The result becomes inaccurate as this quantity grows with 
increased squeezing. Moreover, they did not calculate the squeezed limit of Oijk, assuming it 
to contribute negligibly for an observable bispectrum of C ~ as is the case for near-equilateral 
conhgurations. We give more details on their approach and when it is accurate in §4.2 and 
Appendix D. 

In Ref. [52] Dias et al. took a different approach. They used a next to leading order 
expression in slow-roll for the bispectrum at the time i*, calculated by the same authors in 
Ref. [62]. Such next to leading order expressions automatically provides information on scale 
dependence, and allows for a mild hierarchy of scales. They went on to use their expressions 
for an investigation of the spectral index of the halo bias. The Dias et al. result includes 
the form of aijk in the case of a mild hierarchy, but also takes the form of an expansion 
in I log(A;i/A; 3 )| with only the hrst term calculated, and so again becomes inaccurate as the 
squeezing is increased. Their result is given explicitly in Appendix E. 

A major aim of this paper is to determine when previous expressions can be trusted, 
and when the squeezing becomes large enough that the expressions we present are necessary. 
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3 The squeezed limit of the bispectrum with 6N 

In this paper, therefore, we wish to extend previons work to explicitly calcnlate the bispectrnm 
at time tu in the case of a trnly sqneezed momentnm confignration, ki <C A :2 ~ fes, which will 
involve (perhaps very) different horizon exit times ti <C t 2 ~ ^ 3 - To do so we will employ 
expression (2.2.11), setting T = bnt will fnlly acconnt for the different horizon crossing 
times, and also show how to calcnlate for highly sqneezed confignrations. In contrast to 

the previons section, the new objects we need to calcnlate are Sbq/ci), the field pertnrbation 
two-point fnnction for modes which cross the horizon at time ti evalnated at the later time 
and Q;b(.(/ci, ^ 2 , k^) for the ki ^ k 2 ^ ks confignration. Many anthors have discnssed how to 
propagate the field-space two-point correlation fnnction past horizon crossing [60-62], bnt to 
the best of onr knowledge no one has then employed these techniqnes to explicitly investigate 
the bispectrnm of C beyond cases of mild sqneezing. Moreover, to the best of onr knowledge 

/o'\ 

we are the first to attempt a calcnlation of Oi\-^{ki^k 2 -,k‘i) in the sqneeezed limit nsing a 
backgronnd wave method' which generalizes most easily to a resnlt which is independent of 
a slow-roll assnmption. 

We proceed to discnss how to calcnlate each of these objects now, and then combine 
them to compnte the bispectrnm of ^ in the sqneezed limit for mnltiple field models. 

3.1 Field evolution between different crossing times 

The first task is to calcnlate S^^^(fci). To do this we need to acconnt for the evolntion of the 
field-space pertnrbations on a flat hypersnrface with wavennmber ki between the time this 
wavennmber crosses the horizon and the later time t^. We choose to approach this problem 
in a manner closely connected with the 5N framework, and which therefore provides a nnified 
treatment of the overall problem. It will also make transparent the cases in which analytic 
progress can be made. 

Retnrning to the separate nniverse pictnre, pertnrbations at time are defined by 

64>f\^) = (3.1.1) 

where 0) '(x) is the trne valne, and cj)l ' withont an x argnment is the homogeneons back- 

/o\ 

gronnd valne. In general, for a given i, (/>) '(x) will depend on the valne of all the fields and 
field velocities at the earlier time ti. Assnming, however, slow-roll between ti and the 
fields velocities become fnnctions of the fields. One can write therefore 

(5(/>f^(x) = (3.1.2) 

= + 5</>5-^^(x)) - (l)f\ (3.1.3) 

In analogy with the 5N expression, one can Taylor expand this in the pertnrbation (j(^j-^^(x) 
to give 

^ ,(3) 

+ ... (3.1.4) 

84)) 

'’Although similar results can be found in [12] using a Hamiltonian method and [64] using a second-order 
perturbation theory method. See also [65, 66] where a in the squeezed limit was given in terms of OPE 
coefficients. . 
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In this paper we won’t need the higher order terms. We nse the shorthand 


(3,1) _ d(t>\ 


(3) 


ij ~ 


d(j) 


( 1 ) 


(3.1.5) 


which we will often refer to as the T-matrix’. We wonld like to emphasise that the T-matrix’ 

(T) 

nsed here is constructed nsing the backgronnd cosmology and is analogous to in (2.2.6). 
In Fonrier space we have 


_S- ,(3) T-(3d)xi(l) I 


(3.1.6) 


Objects like the F matrices have been nsed by a nnmber of anthors in the past® [58, 68-74], 
and we note the properties: = 6ij, and so = 6kj for 

arbitrary times ta, th, tc- Moreover since C in (2.2.5) is independent of T, at first order we 
find that 


= ( 3 - 1 - 7 ) 

and so we mnst have 

^(b) ^ Wa)r(a,b)^ (3.1.8) 

X J J L v / 

We can now nse the F matrices to relate correlation fnnctions at different times, and in 
particnlar nsing Eq. (2.3.1) we find 


(3.1. 


('3'\ 

While at one level Eq. (3.1.9) has simply swapped the nnknown correlation matrix {ki) 

(3 1) ^ 

for the nnknown F)^’ ^ matrix, the point is that the problem has been rednced to a matter of 

solving the backgronnd cosmology, which mnst be done anyway to apply the 6N formalism. 

Moreover, in all models where the derivatives of N can be calcnlated analytically, analytic 

expressions shonld also exist for the reqnired F matrix, as we will see in an explicit example 

below. On the other hand, if nnmerical tools are needed to calcnlate the derivatives of N, 

similar tools can be applied to calcnlate The formnlation in terms of F is therefore 

very convenient for explicitly stndying the sqneezed limit. 

3.2 The field-space bispectrum in the squeezed limit 

The next step is to calcnlate o:-jl(ki, k 2 ,k^) in the sqneezed limit. To do so we can adapt an 
idea nsed by Maldacena [1] and snbseqnently other anthors [3, 7] in the single field context. 
The calcnlation relies on the following simplifying assnmption: the long wavelength mode ki 
which exits the horizon at the first exit time, ti, can only affect the mnch shorter wavelength 
modes, which at that time are still deep inside the horizon, throngh its effect on the back¬ 
gronnd cosmology. The ki mode shifts the backgronnd field confignration, and hence changes 
the backgronnd cosmology which is felt by the two modes which exit later - prodncing a 

®r can also be generalised to be fc-dependent, as in [67]. 


= cyr'; 


_ ■p(3,l)-p(3,l)y^(l) /1 \ _ t^{3,1)t^(3,1) 


jk 


*7 


jk 


2kf 
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correlation between the ki mode and the two-point function of the k 2 , k^^ modes in the shifted 
background. 

in the single field case the calculation is performed directly with the modes of Q in the 
comoving gauge, and the shift in background cosmology induced by the long wavelength mode 
exiting the horizon can be thought of as a shift in the time at which the short modes exit. 
This leads to the three-point correlation function of C being related to the tilt of the two-point 
function of the short modes, and to the famous consistency relation of Maldacena®. 

in the multiple field case we are calculating the three-point function of field perturbations 
in the flat gauge, and the shift in the background cosmology induced by the long wavelength 
mode exiting the horizon can be thought of as a shift in the background field values. This 
won’t result in a simple consistency relation for (^, but does allow analytic progress to be 
made in calculating a in the highly squeezed limit. 

To calculate a we follow a similar calculation to that in Cheung et al. [7] for (CCC) iii 
the single field case. We show a derivation of the squeezed limit of a which relies on slow- 
roll - though this derivation can be easily generalised to relax the assumption of slow-roll. 
In the end the final result we get for the squeezed limit of (CCC) must rely on a slow-roll 
approximation at the time of the last horizon crossing, t^, because we will be using the 6N 
formalism at that time. 

We begin our calculation in position space, denoting short-wavelength perturbations 
with a superscript S, and long-wavelength perturbations with a superscript L. We ask how a 
short-wavelength two-point function at t^ is affected by a further long- 

/ o\ 

wavelength perturbation 54>^^ (x). To proceed, we evaluate the long-wavelength fluctuation 

/ o \ 

(x) at the midpoint x = x+ = (x 2 -|- x^)/2. Then we Taylor expand the short- 
wavelength two-point function around the value it would have in the absence of the long- 
wavelength perturbation, denoting this value by a subscript 0. One finds 










(3.2.1) 


... 


where the subscript , m denotes the partial derivative with respect to the background field 
4>m- At this stage we employ the simple soft limit argument discussed above and assume 
that in the squeezed limit the three-point function in momentum space will receive its largest 
contribution from the correlation between the long-wavelength mode, which effectively shifts 
the background cosmology, and the short-wavelength two-point function in the shifted back- 

®Note that the consistency relation holds for all single field models with a Bunch-Davies initial state and 
where the classical solution is a dynamical attractor. 
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ground. This leads to 




v,l{3) 


(x+) 


(3.2.2) 


(xi {x2)S(j)^k'’ (xs)),m ) 

0 


(3.2.3) 


(x+)) { 5 (j)^f^ {yi2)6(J)^k^ (xa)),^ 


0 

(3.2.4) 


[ P ^ ^ ip-(xi-x+)+tq-(x2-X3)v(3)/ ^v(3) / N 

J ( 27 r )3 ( 27 r )3 


0 

(3.2.5) 

In what follows we will drop the subscript |o for notational ease. Now we insert 1 = 
f (i^ki(5(ki + p) to get 




d^ki d^q 


^-.k.-x.-ip.x++i,.(„-x.)(2^)3j,^^ ^ p)sS(p)sfT{<;). 


(3.2.6) 


, (27r)3 (27r)3 (27r)3^ 

Changing the integration variables from p,qto k 2 = ks = + q we find 

(xi) f ^ (X 2 ) <5(/>^i^^X3 ) ) 


d3ki d^k2 d^ka __ 

(27r)3 (27r)3 (27r)3 


iki'Xi—ik2X2—iks'Xa 


( 27 r) 3 d(ki + k 2 + k 3 )sS(fci)S)^,(^(/^ 3 ) 


43 ) 


jk,m^ 


(3.2.7) 


where we have used the relation q = k 3 (l + ©(fei/fea)). From (3.2.7) we can now read off the 
squeezed limit of the momentum space three-point function of the field perturbations 

(2>r)3«(k, + k2 + k3)s2(J:,)S®„(k3) (3.2.8) 


and so 

fciS k2 ^ 2 ’ ^ 3 ) ~ 


(3.2.9) 


This is a very general expression for the squeezed limit of a, independent of the multiple 
field model, relating the squeezed limit of the three-point function of the field perturbations 
to the two-point function of the field perturbations and its derivatives with respect to the 
background fields. It is one of the principal results of this paper. It only relied on slow-roll at 
time fa in (3.2.1) where the Taylor expansion was done only in terms of the field, rather than 
both the field and its velocity. If one promoted the indices i to run over both fields and field 
velocities, then (3.2.9) with this more general index notation would still hold, independent of 
whether slow-roll was valid at time ta. 

At this stage we can make our result more explicit in the case where slow-roll is valid 
between the horizon exit times ti and fa- In this case, we can use the T evolution between 
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crossing times, which by (3.1.9) and (2.3.1) gives 




r(3,i)r(3d)^^.^ 




Hi^r 


ml 


2kf 


2k^ 


^(3,1)^(3,1), HW F(3) 

^ il ^ ml ^jk 


(3) 

m 


(3.2.10) 


2kf 


2k^ 


dN 


where the final line nses the slow-roll eqnations of motion for the backgronnd fields, and for 
clarity we note that dN = d log a is the local measnre of e-folding time (and is not related to 
the 5N formnla). 


3.3 The squeezed limit of the bispectrum of ( 

We now have all the ingredients we need to calcnlate the bispectrnm of C (2.2.11) in the 
sqneezed limit for mnltiple field inflation, and since we now nse the 6N formalism, the resnlts 
which follow in this and snbseqnent sections are only valid in slow-roll models. The schematic 
pictnre of onr approach is shown in Fignre 1. 



p(3,l) 

^ jrn 




Figure 1: Schematic picture of 6N evolution from time T = onwards, using Sbq/ci) from 

Subsection 3.1 and Oi\j^{ki, k 2 , k^) from Subsection 3.2. The F-evolution of field perturbations 
occurs between ti and The blue line is the comoving Hubble radius, the solid red line is 
the inverse of the squeezed wavenumber, fci, while the solid black lines are the inverses of the 
other wavenumbers, /c 2 and k^. 
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Putting Oifjliki-, k 2 , k^) from (3.2.10) and T,'l^\ki) from (3.1.9) into (2.2.11), with T i—)• ta, 

gives 


lim B(;{ki,k2,ks) 

ki<^k2,k^ 


iv^ i\j i\j ^ mk 


/7(3)2 ^(^ 0 ) 

2kl 2^3 dN 




^im ^ jm ^u3 ol3 ^ 



2 " 

(3.3.1) 


and since ki <C ^3 this can be simplified to 



(3.3.2) 


which is one of the main results of this paper. In Appendix A we check that it reduces to the 
Maldacena result [1] in the single field limit. In Appendix B we check that our result for a 
in (3.2.10) agrees with the result of Seery & Lidsey [47] if we take a near-equilateral limit. 
We can then form the reduced bispectrum in the squeezed limit 


, ^h'Liki,k2,k3) 

ki<^k2,k3 5 




+ 




(3.3.3) 


For convenience we define 

(3.3. 
(3.3. 

so that superscript [a] labels the contribution from the term involving a. 

As we noted in §2.3 for the near-equilateral configuration, the term coming from a can 
only be of order slow-roll [47], and so if the reduced bispectrum of C is to be sufficiently 
large to be observable by present or next generation experiments, it will be dominated by the 
second line of Eq. (2.2.11) [63]. 

For our highly squeezed case it’s not immediately clear whether the contribution of 
f^^{ki, k 2 , fca) must be of order e. One might think that because at least one component of 
the vector Aj is of order and because as many F matrices appear in the numerator 

as in the denominator, that the likely order is indeed e. However, while this is likely the 
most common outcome, because is a matrix and appears with different contractions 

in the numerator compared with the demoninator, it is not impossible that the contractions 


(3) 

m 


c \ — aj'( 3) Ar(3) AT(3)-p(3,l)p{3,l) ^ ^ d(j)' 




^ik ^mk 


d(l)m 


l\n k iq i-nq 
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may conspire to make ^ 2 ; ^s) larger than 0{e). In this paper we do not perform a 

general study, but armed with this explicit expression for ^ 2 , ks) we can consider the 

amplitude of this term on a case by case basis. 

Finally, for completeness, we note that in addition to purely scalar correlations, it is 
possible to calculate the three-point functions involving both the scalar curvature perturbation 
and gravitons in the squeezed limit. In Appendix C we use similar techniques to those 
employed above to find the scalar-graviton three-point functions for multiple field models, 
noting that the graviton-only three-point function will be the same as for the single field case, 
as given in [1]. 


4 Scale dependence 


4.1 Spectral index of the halo bias 

An important quantity for large scale structure surveys is the scale-dependent halo bias, 
6b{ki) [45], which is sensitive to how the ratio of the bispectrum to the power spectrum, 
B(^{ki,k 2 ,kz)/Pc^{ki), scales with the squeezed momentum ki, where /ci <C A :2 ~ fes. This is 
captured by the spectral index of the halo bias, nsb = ngq — {ug — 1), where Ugq is the tilt of 
the squeezed limit of the bispectrum with respect to its squeezed momentum ki: 


lim S^(/ci,fc2,/c3) ~ 





(4.1.1) 


with kg some arbitrary scale and roughly constant. Dias et al. [52] investigated nsb in 
multiple field inflation in the case where | log(A:i/A: 3 )| is of order a few, as discussed further 
in Appendix E. Here we would like to explore the highly squeezed case, and find if there are 
significant differences. 

Our results of §3 for the bispectrum in the highly squeezed limit can be applied to 
calculate nsb for large values of | log{ki/k 3 )\. Differentiating log{kiB(^) in (3.3.2) with respect 
to log ki we find 


^sq 


Ar^3)^0)(^D)^{3) g^D)^(3)2) 


2g(i)^(3,i) _ ^D,1) 


r(3,l) 

^mn 


where we have defined 


r (3,1) _ p(3,l)p(3,l) 
'^ij ^ im jm 


dL 


(3,1) 




— _ 

~ dlogki 


V 


( 1 ) 


H(i) 


V 




I jk 


ik 


jk,l 


■p(3,l) _ 
^ iki — 


d4>, 


^p(3,l) 

(1) ik 


(4.1.2) 

(4.1.3) 

(4.1.4) 

(4.1.5) 


and we have used the slow-roll equations to relate cj)i to in (4.1.4). 

The spectral index of the halo bias requires also the tilt of the power spectrum at ki 


ng-l = 


d\og{klPi^{ki)) 


= -2 


^(3)^(3)r(3,i)r(3.i)^(i) 

7y(3) Aji^) T 

^mp 


dlogki 


(4.1.6) 













where 


andiS]^= ’* 


*j 


y(l)2 V(^)' 

This leads to the expression for the spectral index of the halo bias 


(4.1.7) 

(4.1.8) 



jY,(3)jvf )(Arf y® -F 6Ay 74(3)2) 


2,(i)yd)_^0d) 


'^db 

+ 2 

74(3)2) 

^{3)^(3)^(3,l)r(3.1)^{l) 

/\r(3) T 

r(3,l) 

■^mn 


valid for large values of | log(A;i/A: 3 )|, and is another key result of this paper. In Appendix E 
we show that when we take our exit times to be roughly equal (4.1.9) reduces to the same 
form of (17) of [52], given in the appendix as (E.0.2), and in §5.3 we compare approaches for 
a concrete model. 

We note that if k 2 , k^) is negligible, then our expression simplihes to 


nsb = 


ivf)ivf)Arg) (2e(i)Lg^ 


1 ) 


_ p(3d) 


ju-ls) Tyls)/wls) 

JVrn, JVr l\rn 


L 


(3,1) 


+ 2 


^ 7 ( 3 ) 7y (3) T (3,1) 

7 V m i V n ±Jmp 


(4.1.10) 


4.2 Tilts of the reduced bispectrum in the squeezed limit 

In a similar manner to Refs. [51, 75-77], one can study the scaling of the squeezed limit of the 
reduced bispectrum (3.3.3) with respect to ki and fes. In the squeezed limit k 2 ~ k^^ so one 
can parametrize how this depends on scale by differentiating with respect to ki or ^3 leading 
to the tilts 


X _ 
/nl 


91og|/NL(fcl,/C2,fc3)| 
dlog A 


(4.2.1) 


where X = ki or fcs, with the other momenta held hxed in the derivative. Note that = nsb 
since the ki dependence of /nl is captured by the scaling of the ratio B(^{ki,k 2 ,k^)/P(^{ki) 
with ki- We calculate in (F.0.3) of Appendix F. We note that each of these results are 
different to the scale dependence of the equilateral conhguration of Byrnes et al. [51] given 
in Eq. (2.3.5), since we are working only in the squeezed limit. These authors also considered 
the scale dependence of near-equilateral triangles, writing the wavenumbers as ka = C(ak, and 
varying with respect to k, keeping the aa constant. Their near-equilateral result has the same 
form as (2.3.5), with all *’s replaced by the exit time of a pivot scale kp not too different from 
the k's. Their calculation relies on an expansion to hrst order in | log(/ci//c 3 )| and as a result 
is only valid for a small squeezing. We calculate the k tilt using our approach in Eq (F.0.8) 
of Appendix F, where more details on the expansion of [51] can also be found. We emphasise 
again that all of our expressions can be employed for a large hierarchy of scales. 
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5 Employing the T formalism in concrete models 

In this paper we have advocated the use of the T matrices to allow the bispectrum of C to 
be calculated in the highly squeezed limit for multiple held models. These matrices allow 
us to account for the evolution of the inhationary helds between horizon crossing times, and 
help provide compact expressions for the bispectrum and its scale dependence, including the 
contribution from the held-space three-point function. To be of use, however, we must be 

able to calculate the T matrices in concrete settings. 

(3 1) 

In any given model, we could solve for Tl ’ ' numerically. Either by solving the equation 

of motion it satishes [58] from the initial conditions at ti of Tly = 6ij, which follow from 
its dehnition, or by solving the background equations and implying hnite differences as was 
done for the derivatives of N in, for example. Refs. [58, 78, 79]. However, what makes this 
parametrisation particularly useful is that in any model for which the derivatives of N can 
be calculated analytically, the T matrices also admit analytic solutions. This allows us to 
compute the amplitude of the highly squeezed limit of the bispectrum in such models and to 
compare this against previous squeezed limit expressions, as well as against the amplitude of 
the bispectrum for near-equilateral triangles. 

We therefore calculate T analytically in §5.1 for sum-separable potentials. Then we con¬ 
sider the importance of T in specihc settings. In §5.2 we show that in all single-source models, 
where only one held (which need not be the inhaton) contributes towards the curvature per¬ 
turbation, the effect of the T matrices is to cause the reduced bispectrum in the squeezed 
limit to become independent of the squeezed momentum ki , and that because the bispectrum 
scales with ki in exactly the same way as the power spectrum, the spectral index of the 
halo bias will be zero, and so not observable. Then in §5.3 we consider a specihc multiple- 
source model where more interesting results are possible. In particular, we examine a mixed 
curvaton-inhaton model [43, 80-87] allowing for self-interaction terms for the curvaton. For 
this specihc model with our given parameter choices, we hnd that in highly squeezed cases 
relevant for future observations, the bispectrum is suppressed by the T matrices at a level of 
20% when compared to using the existing expressions of Byrnes et al. [51] and Dias et al. [52] 
for this model. In addition, we hnd that the spectral index of the halo bias is enhanced at a 
level of 20% in this model compared with the results that would be obtained with previous 
expressions which assume a mild hierarchy of scales. 

5.1 Calculating T: sum-separable potential 

In the 5N framework, all models which are analytically tractable have a common feature. 
This is that the inhationary potentials are of separable form, either sum or product separable 
[78, 88, 89], or of the generalised sum-separable form of Ref [90]. This is true not only for 
models in which the evolution is tracked during inhation, but also for models in which the 
post inhationary evolution is important such as the curvaton model. In this work we will 
focus on sum-separable potentials and conhrm that we can derive analytic formulae for the 
T matrices, using similar techniques as those originally used for derivatives of N in Ref. [78]. 

We will initially work with the simple case of a two held model: 4>i X write the 
potential as W{(f),x) = U{(p) + V{x)- The slow-roll equations are then 

= 3Hx = -Vx^ 3H'^ = W. (5.1.1) 


17 


Using these slow-roll equations we have 


d(j) 

dx 


X 




(5.1.2) 


The number of e-folds between a flat hypersurface at time ti and another flat hypersurface 
at time is 


AN = 


and using (5.1.2) this gives 


Hdt = / -^dcf) = - / —# 

70(1) 4 > 70(1) U0 


'tl 


(5.1.3) 


AN = - 


/•0(3) TJ y 

/ -#- / — 

70(1) 7^(1) U;^ 


Trdcl)- / 


(5.1.4) 


We also have, by (5.1.2), that 


1 


/■ 0 ( 3 ) 2 

/ I 

/0(i) U^0 J >x 


(5.1.5) 


To determine we need to hnd the following four derivatives of flat hypersurface helds 

(5.1.6) 


/ 00(^) 90(®) 

(3)1) _ I 00(1) 

I 

\ 00(1) 0x(i) 




For flat hyper surfaces, by dehnition, if we vary our position on the initial slice, then AN 
does not alter. This implies that the derivative of AN with respect to held values on the 
initial hat hypersurface satishes: AN = 0 = ANyi). Employing (5.1.4) then leads to two 
independent equations relating the derivatives. Moreover, differentiating (5.1.5) with respect 
to the held values on the initial hat hypersurface yields two further independent equations 
relating the derivatives. Between these four equations we can then solve for each derivative 
and hence determine 


We begin with AN ^(i) = 0 and AN 

x(i) = 0 giving respectively 


^(3) 5^(3) 

y{3) Qyd3) 

[/(I) 

(5.1.7) 

7y(3) Q(j)W 

y^) (97>(1) 

I'll’ 

[/(3) 

t/(3) QyX3) 

E(i) 

(5.1.8) 

)0 

V{3) 

y(l)' 

Next differentiating (5.1.5) with respect to and 

gives respectively 


1 

1 dx^^'^ 

1 

(5.1.9) 

jj(3) d(j)W 

y(3) Q(j)i3) 


1 

1 Sx*'^^ 

1 

(5.1.10) 

^(3) 5^(1) 
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y-U) Q^{1) 

y(l)- 
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Solving, we find 


p{3,l) 


a<^(3) 90(3) \ 
90(1) 9x(i) 1 
9x(3) 9 x(3) I 
90(1) 9x(i) / 


/ (^(l)+y(3)) 

^(1) iy(3) 

V® ([;(1)_[7(3)) 

1 [/’(I) W(3) 

\ ,0 


(y(i)-v(3)) \ 

T/(l) 

y® (y(i)+i7(3)) 
yW vy(3) / 


(5.1.11) 


which is the analytic calcnlation of T for two-held snm-separable potentials we reqnired. Note 
that i labels the rows and j labels the colnmns. 

For a model with n helds one can calcnlate the F matrix by generalizing the two-held 
case. In what follows we will snspend the snmmation convention of repeated indices, and 
where a snm shonld be taken we explicitly state this. We take the snm-separable potential 

We have, in analogy with (5.1.4), 


AN = - 



w:‘^‘ 


(5.1.12) 


where VF/ = dWi/d(j)i for which we can write the eqnations AiV (i) = 0 in analogy to (5.1.7) 

^<Pj 

and (5.1.8). We also have the relations, in analogy to (5.1.5), for all l,i 



(5.1.13) 


for which we can take derivatives with respect to giving eqnations analogons to (5.1.9) 
and (5.1.10). Combining these eqnations with those from AN (i) = 0, after some algebra we 

arrive at 


(3,1) _ d(j)[ 


(3) 


ij 


d(J)' 


( 1 ) 


pp'(3) (w)^) - 


which rednces to (5.1.11) in the two-held case. 


(5.1.14) 


5.2 Single-source models 

In what follows we will assnme for now that in the highly sqneezed limit we can neglect 
k 2 , ks), the contribntion to the bispectrnm from intrinsic held-space three-point fnnc- 
tion. This needs to be checked on a case by case basis. 

In the single-sonrce case where only one held, which we denote y, contribntes to C, the 
bispectrnm (3.3.2) (withont B^^\ki,k 2 ,ks)) simplihes to 


lim Bi;{ki,k 2 ,ks) 

ki<^k2,k3 


X XX X X'm X'm 


if(3)2 

2kl 2kl 


(5.2.1) 
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Now we can use (3.1.8) to relate to through the T’s to give the bispectrum and 
reduced bispectrum as 


lim B(^{ki,k 2 ,h) 
ki<^k2,k:^ 


lim -/NL(fcl,/C2,/C3) 
ki<^k2,k^ 5 


2wC)iV(3)jvC):f^l!^ 

z,iv^ m 973 ^13 


N. 


( 3 ) 


XX 


2Pc(fei)Pc(fc3)- ,^,2 
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N. 


( 3 ) 


N. 


XX 

( 3 ) 2 ' 


X 


(5.2.2) 

(5.2.3) 

(5.2.4) 


We note that the effect of the T matrices has been to make the bispectrum proportional to 
2 P(^(ki)P(^{k 3 ), so that the reduced bispectrum is independent of ki, and coincidentally of the 
same form as would be derived in the near-equilateral regime assuming t* ~ ts. Because the 
bispectrum scales with ki in exactly the same way as the power spectrum, the spectral index 
of the halo bias will be zero, and so will not be observable unless the intrinsic contribution 
we have neglected is important. 

Note that the results of this subsection did not rely on assuming a sum-separable poten¬ 
tial - all single source models satisfy (5.2.4). In the next subsection we will look at a specific 
multiple-source model which does rely on the assumption of a sum-separable potential. 


5.3 Multiple-source models: the mixed curvaton-inflaton model 

We now consider multiple-source models, and will consider the concrete example of a curvaton 
model for which both the inflaton, ip, and curvaton, y, contribute towards (p. In order to make 
analytic progress we will need expressions for the derivatives of N in this multiple-source 
model to combine with our analytic expression for F, valid for sum-separable potentials. 

More detailed studies of the curvaton scenario can be found in the literature, see e.g. 
[80-85]. Here we just give a very brief account outlining our parameter choices and quote the 
results for the N derivatives. Our main focus will instead be on the differences between our 
new expressions compared to existing formulae in the literature. We will show these diffences 
graphically, where the plots are produced using the expressions for the N derivatives that 
follow, together with the F matrix for sum-separable potentials. 

For the inflaton and curvaton, we take the potential 

W{(p,x) = + ^mlx^ + Xx"", (5.3.1) 

where the curvaton is given a self-interaction with coupling A, and n > 2. If we are to see a 
difference between the bispectrum in the highly squeezed limit and the near-equilateral limit, 
it is natural to expect that we will need the field configuration to differ significantly between 

different exit times, otherwise the dcpi would not evolve between ti and and we would just 
(3 1 ) 

get FF’ ' = 5ij. This means that we expect to see interesting effects for models which generate 
a significantly scale-dependent non-Gaussianity. In the present model, this means that the 
curvaton field, while light, cannot be completely frozen. To achieve this we take the mass of 
the curvaton and inflaton to be the same rriy- = m^. 

The self-interacting curvaton model can give significant scale dependence of /nlj [43, 91], 
and analytic expressions for the derivatives of N are available in the limit of weak self¬ 
interaction Sk = 2Ax^~^/m^ <C 1 [43]. Here, and in what follows, a subscript k indicates the 
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result is a function of the value of the curvaton and/or inflaton helds at the time, when a 
given /c-mode exits. 

The relevant N derivatives are given by [43, 50, 53] 


N^ = 


1 


^ 3 a, 


N - — 


cr, 


osc 

// 


(Tr. 


where aosdxk) oc xa, (^1 + 


+ 

^OSC \ ^osc 

-l/(n-2) 


(5.3.2) 

(5.3.3) 

(5.3.4) 

(5.3.5) 


The parameter rjec denotes the value of 3py./{3p^ + 4:p.y) at the time of curvaton decay, where 
p is the energy density and the subscript labels the species, with 7 denoting radiation. In 
the example which follows we take to be rjec = 0.02. We have neglected which are 

much smaller than for the parameter choices considered. 

We take the initial condition^^ tpQ = 16 which leads to 63.5 e-folds of inflation and assume 
that all scales which exit the horizon after this time are within the horizon today and poten¬ 
tially observable. We will also take xo = 2 x 10“^ to generate a signihcant non-Gaussianity. 
To see a signihcantly scale-dependent non-Gaussianity we take the self-interaction with n = 6 
and A = 0 . 2 , which gives Sk ~ 0.07. 

The power spectrum for this model is shown in Figure 2a, and the tilt of the power 
spectrum is shown in Figure 2b. In order to see a large scale dependence of /nl for this 
simple model, we took the curvaton mass to be as large as the inflaton mass^^. However this 
has the adverse effect of making the tilt of the power spectrum to not be as red (negative) as 
observations suggest. This could be hxed by beginning with a more complicated model, and 
not requiring the curvaton mass to be as large, but this is not our focus here, and instead we 
focus on the results for the bispectrum. 

We can now compare our new expressions with F’s against existing expressions in the 
literature. In the 3d plots, the near-equilateral limit log(fci/A; 3 ) ~ 0 is on the far right of each 
plot, and the squeezing increases as you move away from this corner. The highly squeezed 
limit, where log(A:i/A: 3 ) ~ —20 is on the far left of each plot. 

In Figure 2c we plot the bispectrum /3 = Akfk^Bi^ to show how the bispectrum scales 
with log(A:i/fco) and log(/c 3 //co) where ko = at the initial time to, where we have set 
ao = 1. The blue surface uses our new expression (3.3.2), the red surface uses (96) of Byrnes 
et al. [51] given in (D.0.1), and the yellow surface uses (14) of Dias et al. [52], given in 
(E.0.1). In the highly squeezed limit we see a percentage difference of about 20%. The near- 
equilateral results are within a few percent of each other for each of the three expressions - 
the small difference is because previous authors have included an additional constant term 
which appears at next order in slow-roll in their expression for S at horizon crossing discussed 
in Footnote 5 and Appendix D. 

^°Note we are still working in units where the reduced Planck mass is set to one. 

A potential issue for the mixed curvaton-inflaton model with equal masses Is that since both helds begin 
oscillating at exactly the same time, the curvaton energy density remains subdominant for a long time - see 
[92] for a discussion on this. 
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(c) Bispectrum (3 = in (3.3.2) (d) ^[“1 = from first line of (3.3.2) 



(e) Reduced bispectrum /nl, (2.2.12) 


(f) Spectral index of the halo bias, (4.1.10) 


Figure 2: In the 3d plots, the blue surfaces are our expressions with F’s. Red surfaces are 
Byrnes et al. [51] expressions. Yellow surfaces are the Dias et al. [52] expressions. The 
near-equilateral regime is on the far right of each plot. The highly squeezed limit is on the 
far left of each plot. The green surface shows the [a] contribution is negligible for this model. 


In Figure 2d we plot the contribution to the bispectrum from the held-space three-point 
function /?[“] = where is the hrst line of (3.3.2). We see that this is a factor 

of 10^ smaller than the total bispectrum for all scales considered, and so all the [a] terms in 
observables can be neglected for this model. 

In Figure 2e we plot the reduced bispectrum, where the blue surface is our new expression 
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(3.3.3) and the red surface is (96) of [51] divided by the factor 2P^(A:i)P^(/c3), i.e. the curly 
brackets of (D.0.1). Dias et al. don’t give an explicit expression for the reduced bispectrum, 
so we don’t plot this here. Again, in the highly squeezed limit we see a percentage difference 
of about 20%, and close agreement in the near-equilateral conhguration. 

Finally, in Figure 2f we plot the spectral index of the halo bias. Blue is our new expression 
(4.1.10), and the yellow surface is (17) of [52], given in (E.0.2). Byrnes et al. [51] don’t give an 
expression for the spectral index of the halo bias so we don’t plot this here. Similarly to the 
other observables we see a percentage difference of about 20% in the highly squeezed limit. 
We see agreement to within a few percent in the near-equilateral limit, the small discrepancy 
being due to the Dias et al. expression being evaluated at the exit time of kt ^ ki + k 2 + k^ 
compared to ours being evaluated at t^. 

The 20% level difference for all these observables for a squeezing of | log(A:i/A: 3 )| ~ 0(20) 
can be estimated heuristically as arising from a scale-dependence of the reduced bispectrum 
of 0(0.01), multiplied by the squeezing. We thus expect similar levels of discrepancy for any 
model where the scale-dependence of the reduced bispectrum is of similar order to the scale- 
dependence of the power spectrum, when considering a squeezing of | log(/ci/fe 3 )| ~ 0(20), and 
that the percentage difference will scale linearly with the scale-dependence of the bispectrum. 

The expected levels of squeezing for three future experiments were shown in Table 1 of the 
introduction. The results of this section show that inclusion of the effects of held evolution can 
be important when computing the theoretical predictions of a model for comparison against 
observations for large squeezing, even for this simple model. 

6 Conclusion 

In this work we calculated the squeezed limit of the bispectrum of the curvature perturbation 
for multiple held inhation. Different scales involved in one triangle of the bispectrum will 
exit the horizon at different times, and previous analytic expressions have been limited to a 
mild squeezing where the exit times are roughly equal. Observations can at present probe 
only a mildly squeezed limit, but future large-scale structure surveys and observations of 
CMB //-distortions will be able to probe a highly squeezed limit. It is important, therefore, 
to have accurate theoretical predictions for this highly squeezed limit in order to ensure the 
uncertainty in the prediction is less than the uncertainty in the data. For certain models, 
our results give a correction at a level of 20% in the highly squeezed limit compared to 
extrapolating existing expressions, valid only in the mildly squeezed limit, to the highly 
squeezed limit. 

In order to study this highly squeezed limit, we suggested using the elegant F matrix 
formalism to account for the evolution on superhorizon scales of the perturbations between 
exit times. We also calculated the intrinsic three-point function of the held perturbations, 
a, in the highly squeezed limit for the hrst time. We did so by appealing to a soft limit 
argument, previously used in the single held context for the curvature perturbation. Together 
these elements allowed us to extend 5N expressions for the bispectrum of C to account for 
multiple crossing times. From this expression, we then obtained the reduced bispectrum 
and the spectral index of the halo bias. Working with a specihc model, the mixed inhaton- 
curvaton scenario with a self-interacting curvaton, we checked the difference in our theoretical 
prediction, valid for a large squeezing, against existing predictions, valid for a mild squeezing. 
As would be expected we found signihcant differences, especially in the highly squeezed limit 
for the cases in which there is signihcant scale dependence in the reduced bispectrum. 
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The overall aim of this paper was to provide clarity in how to confront models of inflation 
against observations sensitive to the squeezed limit of the bispectrum. Our results however, 
could also be useful to check numerical methods in the squeezed limit. From a theoretical 
and observational point of view, soft-limits - of which the squeezed limit is the simplest 
example - are of considerable interest, and in future work we hope to consider soft limits of 
higher n-point correlation functions for multiple held inhation, using a similar approach. For 
n > 3, the story can be more interesting than for the bispectrum, where the only soft limit 
is when a single external momenta becomes small. Firstly, one can consider multiple-soft 
limits, where more than one momentum becomes smaller than the others. Moreover, one can 
also consider the collapsed limit, when an internal momentum becomes soft. Observing these 
higher-point functions may be even harder than for the bispectrum, but it is still important to 
have theoretical predictions for multiple held inhation to constrain models using observational 
limits, in particular to search for deviations from single held inhation. 

Finally, we mention that the intrinsic term in the bispectrum coming from a was neg¬ 
ligible compared to the other term in the case study presented. Calculating this expression 
explicitly allowed us to determine this, but more work is required to investigate whether this 
term can ever be as large as, or even dominate over, the other contribution. 
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A Reduction to single field case 


We should check that our expression (3.3.2) for reduces to the Maldacena single held 
squeezed limit result [1] in the case of one slowly rolling scalar held, given by 


lim BQ{ki,k 2 , ks) ^ -(nf^ - l)P(;{ki)P(^{k 3 ) (A.0.1) 

ki<^k2,k3 


with Ug — 1 = 2r]v — 6ey, and the superscript (3) denoting evaluation at time when k^ 
exits. Here ey = \{V'/V)‘^ and rjy = V"jV are the potential slow-roll parameters. 

Thus we begin with a slowly rolling single held (j). In this case we have 




H _ 1 

ij) \/2ey ’ 




2 ey - rjv , 


d Id 
d(j)^^'> (^(3) dts 


(A.0.2) 


and the T matrix is just the number /Nj^\ Then (3.3.2) 


simplihes to 
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as required. 


B Recovering Seery &: Lidsey result in near-equilateral limit 

A further important check is that the near-equilateral limit of our result for a, (3.2.10) goes 
over to the result of Seery & Lidsey [47], given earlier in (2.3.3), valid when ki is small, but 
not so small as to change the exit times appreciably. That is, we want to check that 


livD. = lim 


tl—>-£3 


27r 


ki<s:k2,k3 kfk^k^ 

^ 4Lr(3 )ITT (^1 + 2^3) + 7,h - kik2 
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6 perms 
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(B.0.1) 


Beginning with the RHS, one can do the sum over all six permutations, then take the 
slightly squeezed limit to get the RHS equal to 




(3) 


^ ^ Akfkl dN 


(B.0.2) 


/Q'\ 7o l'^ 

which is exactly the ti —)• ts limit of in (3.2.10), in which Tb’ —?• 5ij, contracted with 




C Squeezed limits of graviton correlators 

Three-point functions involving gravitons (tensors) are likely significantly harder to detect 
observationally than those just involving scalars. Nonetheless they are interesting to calculate 
with a view to observations in the more distant future, and from a theoretical perspective. 

Maldacena found squeezed limits of scalar-graviton and graviton-graviton three-point 
functions in the case of a single scalar field [1]^^. Here we use our soft-limit argument to 
calculate these in the multiple field case. 

First order graviton perturbations, denoted by 7 , are gauge invariant in contrast to 
scalar perturbations. They are defined as the transverse traceless perturbations of the spatial 
metric, /i/j, such that 


hij = + 2C)6ij+ 'yjj] (C.0.1) 

where 7 // = 0 = di'jjj. We can Fourier expand 7 as 

/ d^k 

where s, r indices label the polarization of the graviton, and the polarization tensors efj(fc) 
satisfy ejj{k) = 0 = kiejj{k) and ejj{k)e}j{k) = 26sr- The two-point function of the graviton 

^^Maldacena and Pimentel [93] found graviton-graviton correlators for gravity theories not restricted to 
Einstein gravity, using a de Sitter approximation. Squeezed limits of correlation functions involving gravitons 
in models of quasi single field inflation were considered in [94] 
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is given by 


(7j.7fa) = (2»)=’«(ki+k2)Fyl‘>(»:i) 


= S. 


2 H^^y 


(C.0.3) 

(C.0.4) 


We now consider the squeezed limit of three-point correlation functions involving gravi¬ 
tons. When the soft mode ki is that of a C) there will be a correlation between two 

7 ’s by way of a similar soft limit argument applied now to a ^77 correlator, giving the result 

<=.S,fa*^'“^k,7E3> = (27)=‘<5(ki + k 2 + k3)JvP>rjy>r»;‘>s7(ti)77>(ks) (c.o.5) 


which can be contrasted with the corresponding single field result given in [1], In the single 
field case, the result is proportional to the tilt of the graviton power spectrum, providing 
another consistency relation between observables. Now in the multiple field case, this consis¬ 
tency relation no longer holds, but instead the squeezed limit three-point function is related 
to how the two-point 7 correlator depends on the background scalar fields (j)i. 

When the soft mode is instead a graviton, we can refer to Maldacena’s argument that 
when the Ck 2 )Ck 3 modes exit, the graviton with momentum ki exited much earlier and is 
already frozen, so that fluctuations of C at time ts will be those in the deformed geometry of 
the background 7 kmode. The main effect of the deformation of the background geometry 
is to change the djjk^k^ —>■ Sjjk^k^ — Jijk^k^ inside the correlation function of the two C’s 
(equivalently in the second order action for Cfca)- Putting this into a soft limit type argument 
gives 

hm^jy^^Ck^Cka) « -(2vr)3<i(ki + + ks)P;''^^\ki)e}j{ki)kiki (C.0.6) 

Similarly, the squeezed limit of the three-point correlator is exactly as given by Maldacena [1] 

h«te.t,<^k.7E37E.)»-(27fi(ki+k2 + k3)P7»(ti)e»„(t,)k't3^^p;‘'3l(t3). (C,0.7) 

We highlight that the result (C.0.5) may lead to new shape dependence compared to the 
single field results. 


D Reduction to Byrnes et al. 


It is also important to check our expressions match those which have previously appeared in 
the literature in the near-equilateral, midly-squeezed configuration, which is the overlapping 
regime of validity. A result for the squeezed limit of the bispectrum was given by Byrnes et 
al. in Eq. (96) of Ref [51] 
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where c = 2 — log2 — 7 , with 7 the Euler-Masheroni constant, and Uij is given by (4.1.8), 
and since mj is of order slow-roll, it can be evaluated at any time. Their result is valid 
for small log{k 3 /ki)uij with the intrinsic contribution from the three-point function of held 
perturbations neglected. It was derived by considering the two-point correlation at unequal 
times, calculated using the expression 


= 


6 ij -h log 

Ki 


Hj 




(D.0.2) 


which is valid for small log(/i; 3 //i:i)ttjj. However, for a large squeezing, \og{k‘i/ki)uij may not 
be small, even though Uij is of order slow-roll, and one will instead need the full T expression 
(3.1.6) for the evolution of the held perturbations. The T matrix can be written formally as 
a time-ordered exponential [58, 68-70, 72, 73] 


p(3,l) 

b 


T exp 



Uij{N)dN 


(D.0.3) 


where A^i is the number of e-folds corresponding to time ti and similarly for N 3 . Note that 
in the limit of small log^k^/ki)uij, we have, at leading order in log{k 3 /ki)uij 


rSf).sh,,+log(|)4') + ... (D.0.4) 

As discussed in Footnote 5, these authors used a next-order in slow-roll expression for S, 

^ + 2™b)- (D.0.5) 

Taking our expression for the squeezed limit of the bispectrum (3.3.1), and substituting the 
RHS of (D.0.4) for T, and replacing 5ij/2ki with the RHS of (D.0.5), we recover (D.0.1). 

The reason we didn’t need to use (D.0.5) in the main part of this work was because we 
were throughout working to leading order in slow-roll, rather than next to leading order in 
slow-roll. The only time we needed to consider slow-roll terms, such as Uij, are when they ap¬ 
pear multiplied by log(A; 3 /A:i), which can be as large as 0 ( 20 ), in which case | log{k 3 /ki)uij\ ~ 
1. Note that the expansion in Eq (D.0.4) cannot be truncated for \ \og{k 3 /ki)uij\ ~ 1. In 
the highly squeezed limit this is why the full expression for T given in Eq (3.1.5) needs to be 
used instead of Eq (D.0.4), even though we can safely neglect the slow-roll correction to the 
power spectrum in Eq (D.0.5). 


E Reduction to Dias et al. 


Dias et al. [52] used a next-to-leading order expression for the bispectrum [62] 


lim B(^{ki,k 2 ,k‘i) 

ki<^k2,k3 
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to calculate the spectral index of the halo bias as 

^ dlogki ^ ^ NiNi 

(E.0.2) 
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where the right hand sides of both expressions are evaluated at the time tt when fc* = fci + 
k 2 + ^3 crosses the horizon, and where 
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in which all quantities on the right hand side are again evaluated at tf, and 7 is the Euler- 
Masheroni constant. These results can be trusted for a mild hierarchy of scales, where 
I log(/ci//c 3 )| is of order a few. We now check that our expression for nsb, (4.1.9), can re¬ 
cover the Dias et al. result, (E.0.2), in the limit where the exit times are very close. To do so 
we begin with the expressions 
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Now assuming we can swap the limit of differentiation with respect to ki and the limit of 
taking ti —)> and using (D.0.4) we get 
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Substituting this into nsb of (4.1.9) and setting T^y ^ —)> 6ij we get 
^(3)^D)(^(3)y(3) 

lim r? r _ _ 9 ^^ \ h ,j _gj_ / ij ^ i _j_ lj 
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which is the of the same form as (E.0.2) when (E.0.3) and (E.0.4) are substituted in. Note that 
in their expression everything on the RHS is instead evaluated at exit time of kt = A:i-|-fc 2 + fc 3 , 
rather than ta, but in the limit where the exit times are very close, this won’t affect the result 
signihcantly, and we recover their result. 


F Tilt of reduced bispectrum in the squeezed limit 


As discussed in §4.2 one can study the tilts of the reduced bispectrum, /nl, in the squeezed 
conhguration with respect to any combination of the /c-modes which it involves. In particular, 
one can calculate how /nl of (3.3.3) varies with respect to ki ~ ^2 or k^,, or some combination 
of them. The dependence can be parametrized by 


X = c^IosI/nlI 
d log X 
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for X = fci,/c 3 . In §4.2 we found = nsb, where ngb was calculated in §4.1 in (4.1.9). 
To find we write (3.3.3) in a form where the second square bracket contains all the 
dependence 
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Note that we have neglected the intrinsic contribution in Eq. (3.3.3) for simplicity. 

To compare with observations, one might wish to use the variables in [95], given by 
k, a, 13, defined as 
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which in the squeezed limit are related to ki , k^ by 
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with d ~ 0. We can use the chain rule to calculate 
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Note that in the squeezed limit /3 ~ 1, and so if we use our expression for the squeezed limit 
of /nl in (3.3.3), we shouldn’t vary jS significantly away from 1. 
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